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Abstract 

Even though the solar wind is highly supersonic, intense ultra-low frequency (ULF) wave activity 
has been detected in regions just upstream of the bow shocks of magnetized planets. This feature 
was first observed ahead of the Earth's bow shock, and the corresponding region was called the 
ULF wave foreshock, which is embedded within the planet's foreshock. The properties as well as 
the spatial distribution of ULF waves within the Earth's foreshock have been extensively studied 
over the last three decades and have been explained as a result of plasma instabilities triggered by 
solar wind ions backstreaming from the bow shock. Since July 2004, the Cassini spacecraft has 
characterized the Saturnian plasma environment including its upstream region. Since Cassini's 
Saturn orbit insertion (SOI) in June 2004 through August 2005, we conducted a detailed survey 
and analysis of observations made by the Vector Helium Magnetometer (VHM). The purpose 
of the present study is to characterize the properties of waves observed in Saturn's ULF wave 
foreshock and identify its boundary using single spacecraft techniques. The amplitude of these 
waves is usually comparable to the mean magnetic field intensity, while their frequencies in the 
spacecraft frame yields two clearly differentiated types of waves: one with frequencies below the 
local proton cyclotron frequency (Qu+) and another with frequencies above Qh+- All the wave 
crossings described here, clearly show that these waves are associated to Saturn's foreshock. In 
particular, the presence of waves is associated with the change in 6 Bn to quasi-parallel geometries. 
Our results show the existence of a clear boundary for Saturn's ULF wave foreshock, compatible 
with 6 Bn ~ 45° surfaces. 
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1. Introduction 

When the supersonic solar wind plasma from the Sun encounters an obstacle, a bow shock 
is formed. The incoming solar wind particles (electrons and ions) upstream from the bow shock 
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have no information about the obstacle, except in a region magnetically connected to it. This re- 
gion is known as the foreshock. At the bow shock, a small fraction of the solar wind particles are 
reflected in the sunward direction. These backstreaming particles are subjected to the solar wind's 
ExB drift, where E= -v sw xB/c, is the solar wind's convective electric field, B is the interplan- 
etary magnetic field (IMF), v sw is the solar wind velocity and c is the speed of light. The ExB 
drift velocity is the same for all backstreaming particles, and perpendicular to the IMF. As a result, 
the guiding centers of all backstreaming particles move within the v sw -B plane, gradually drifting 
away from the field line tangent to the bow shock toward the inner part of the foreshock and being 
segregated according to their parallel velocities. Electrons, because of their much smaller inertia, 
are much less affected by this drift and their presence can be detected right next to the field line 
tangent to the bow shock. The combination of solar wind and energetic bac kstreaming electrons 



results in the producti on of Langmuir waves at the electron plasma frequency (|Fuselier et al.l . ll985 



f i ink 

Sigsbee et al., 2004a.b). Backstream ing ions on the other hand, can drive a number of plasma in- 



stabilities (|Garyl . 1 19931 : IConvery and Gary . 1997 ). leading to the generation of waves. The ion 
foreshock is then characterized not just by the presence of a small fraction of backstreaming ions, 
but also by the generation and propagation of plasma waves around the local ion cyclotron fre- 
quency. 

The understanding of planetary foreshocks is far from complete. The most studied case is the 
Earth's fore shock, which is reasonab l y wel l understood thanks to single and multi- spacecraft mea- 
surements (|Tsurutani and RodriguezLll981l : lEastwood et al. . 12005b. The first observations from the 
Earth's foreshock were made by the dual spacecraft ISEE (lOgilvie et all Il977|) which identified 
different types of backstrea ming ion distributions: reflected (now cal led field-aligned beams), in- 
termediate, and diffuse (|Gosling et al.l . ll978l : IPaschmann et all 1 1981). The classification of back- 
streaming ion populations into these three types was made on the basis of two-dimensional velocity 
distribution functions and energy-time spectrograms. Further results from ISEE demonstrated the 
existence of gyropha se-bunched and gyrotropic backstreaming ion distributions in the foreshock 
dGurgiolo et all 1 1983b . 

The field-aligned distributions are ty pically observed at and near the leading edge of the ion 
foreshock without the presence of waves (|Paschmann et al. . 1979 ). Behind the field-aligned beam 
region, gyrophase-bunched distributions are detected, while diffuse distributions are found even 
farther away (i.e. downstream) from the ion foreshock boundary. The asso ciation of different 
linear and non-linear waves to different ion distributions was first studied by Hoppe and Russell 



19831) . Gyrophase-bunched and diffuse distributions are in fact observed in the presence of ultra- 
low frequency (ULF) waves. In particular, gyrophase-bunched distr ibutions coexi s t with ULF 
quasi-monochromatic waves with substantial amplitudes (SB/B « 1) (IMazelle et all 120031) . The 
production of gyrophase-bun ched ions and associated ULF wave s has been studied numerically 
( Hoshino and Terasawa . 1985). theoreti c ally (Mazelle et al. . 2000l) and observationally, using data 



from WIND spacecraft (Mezian e et all |2001|) . On the other ha nd, non-linear, steepened waves 



have been found to be associated with diffuse ion distributions dlloppe et al. . 1981 ). 

Evidence of foresh ocks has also been found in other planets as well. In the case of Mars, 
Trotignon et al. I d 1992b found electrostatic waves which are polarized along the interplanetary 



magnetic field, and the ir peak intensity occurs at or near the local solar wind plasma frequency. 
Greenstadt et al. I dl987b identified the ULF wave foreshock boundary of Venus through obser- 



vations made by the PVO (Pioneer Venus Orbiter) magnetometer. Plasma and magnetic field 
observations from the Voyager 2 spacecraft reveal ULF waves in the solar wind, which are asso 



dated with Neptune's, Jupiter's and Uranus' foresho cks (see for example Belcher et all 11991 



Bavassano-Cattaneo et all 1 19871 : iRussell et all 11990 ). Plasma waves have also been detected 



upstr eam from Mercury's bow shock using Mariner 10 me asurements ( Fairfiel d and Behannonl . 



1976|) . A brief review on Mercury's foreshock was made by lBlomberg et al.l (|2007|) . 



The region of ULF wave activity is embedded in the ion foreshoc k, and its boundary is known 



as the ULF wave foreshock boundary. In the case of the Earth, iDiodato et al.l 41976) located 
this boundary in a statisti c al way by using combined magnetic field and plasma data from Heos 
1. Greenstadt and Baum (1986), using ISSE 1 data introduced the so called solar foreshock 
coordinates, and determined the position of the ULF wave foreshock boundary. For different 
IMF cone angles (they split their study into 6 Bx = 20° - 30° and 9 Bx = 40° - 50° data sets), 



Greenstadt and Bauml (LL986) obtained different orientations of the ULF wave foreshock bound- 
ary. Other studies have shown that for 6 Bx ^ 45° cone angles, th ere is a well defined u pstream 
boundary, and this boundary intersects the bow s hock at 6r„ » 50° (ILe and Russelllll992a|) . Using 
multi-spacecraft data obtained by Cluster, Eastwood et al.1 (120051) presented two case studies on 
the onset of ULF foreshock waves, relating their appearance to changes in the orientation of the 
IMF. The increased wave activity is associated with the change in 6 Bn to quasi-parallel geometries 
(0 Bn < 45°). These findings are consistent with previous single spacecraft studies. 

Although the observations at Earth represent a vital element in the study of the physical pro- 
cesses occurring on planetary foreshocks at large, these phenomena necessarily occur at the partic- 
ular parameter values relevant for our planet. To explore how these processes change in parameter 
space, it is just as important to make in situ observations around other planets. For instance, the 
solar wind properties vary with heliocentric distance. The Parker angle between the stream lines 
and the radial directio n to the Sun at 1 AU is pr edicted to be about 45° for a constant solar wind 
velocity of 429 km/s ( Thomas and Smith , 1980h. At Satu rn's distance (approx. 9 AU), the Parker 
angle is predicted to be about 85°. Jackman et al. |2008) examined in detail the hourly averaged 
IMF data provided by Cassini from 13 August 2003 to 14 November 2004, and found an average 
Parker angle equal to 86.8°+0.3° for a solar wind speed of 500 km/s. 

The location and shape of planetary bow shocks are determined by the properties of the solar 
wind flow and by the size and shap e of the obstacle to the flow. Saturn's bow shock was first 
observed by Pioneer 11 in 1979 (|Acuna and Nessl . Il980|) . The firs t crossings made by C assini 
during Saturn Orbit Insertion (SOI) were analyzed and discussed by lAchilleos et al.l (|2006|) . They 
presented evidence of magneto spheric compression during Cassini's first immersion into the mag- 
netosphere and the properties of the solar wind upstream from Saturn's bow shock for the first 
six bow shock crossings observed by Cassini. The magnetic signatures of these bow shock cross- 
ings showed a clearly defined overshoot and foot regions associated with the quasi-perpendicular 
geometry (6 Bn > 45°). Usin g magnetic field and plasma observations made by Cassini between 
June 2004 and August 2005. Masters et al.l (|2008|) presented a static model of Saturn's bow shock. 
The model was obtained by fitting a conic section to the first 206 crossing s observed by Cas sini. 
Based on observations from Pioneer 11, Voyager 1, Voyager 2 and Cassini IWent et al.1 (|201 1 ) de- 
rive a small eccentricity for Saturn's bow shock and found variations in the shock subsolar distance 
associated to variations in the solar wind dynamic pressure. 
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Using the Magnetospheric Imaging Instrument (MIMI) on board Cassini. lKrimigis et al.l(|2005l) 
showed measurements of the energetic ion population upstream from both the dusk and dawn sides 
of the Kronian magneto sphere. During the approach phase and first orbits of Cassini, the obser- 
vations revealed the presence of a series of distinct upstream bursts of energetic hydrogen and 
oxygen ions up to distances of 120 Saturn's radii. They concluded that these oxygen upstream 
events m ust be particles leakin g from Saturn's magnetosphere under favorable IMF conditions. 
However. iThomsen et al.1 (I2007|) studied 45 hours of mass-resolved observations of Cassini Plasma 
Spectrometer (CAPS), which were performed upstream from Saturn's bow shock. The observa- 
tions show supra-thermal ions composed of H + and ions with m/q = 2, presumably solar wind 
He + + , with no detectable cont r ibutio n from magnetospheric water group ions. 



Bavassano-Cattaneo et al.1 (| 199 II) reported the first evidence of upstream low-frequency waves 



in Saturn. During these observations, the spacecraft was magnetically connected to Saturn's bow 
shock. Their results suggest that these waves were associated to the planet's foreshock. These 
waves displayed a period of 550 seconds in the spacecraft frame and a relative amplitude of 0.3. 
Also in the spacecraft frame, the waves are left and right-hand elliptically polarized, and propagate 
at about 30° with respect t o the ambient magnetic field. During the first three orbits of Cassini 
spacecraft, Bertucci et al.1 (120071) present a characterizati on of low-frequency w aves associated 
with Saturn's foreshock based on Cassini magnetometer (Dougherty et all l2004f) . As a result of 
their survey, they identified two distinct types of waves. They found a large majority of waves 
with spacecraft-frame frequencies below the local proton cyclotron frequency (f2n+ = eB /m H c). 
These waves are ph ase-steepene d and display a left-hand elliptical polarization as seen by the 
spacecraft. iBertucci et al. (2007) int erpreted these waves as fa s t ma gneto sonic waves. This kind 
of waves is the same presented by iBavassano-Cattaneo et al.l (|199l|) . In a second group, they 
found waves with frequencies above Q H +, quasi-monochromatic and steepened with a right-hand 
circular polari z ation, propagating at small angles with respect to the ambient field. According to 



Eastwood et al. 



Bertucci et al. (2007) these waves could be Alfven waves similar to those observed at Earth by 



(120031) . 

In this paper, we have used Cassini magnetometer data and single-spacecraft techniques to 
study the morphology of Saturn's ULF wave foreshock. In situ observations made by the Vector 
Helium Magnetometer (VHM) on board Cassini are presented in section |2] A brief study of ULF 
waves in the magnetic field associated with Saturn's foreshock is presented in section|3j In section 
H] we show, for the first time, the determination of the outer boundary of Saturn's ULF wave 
foreshock. For this purpose we identified a large number of crossings of Cassini to (or from) the 
wave region. The selection criterion for these crossing is described in section 14. II In secti ons I4~2l 
and [43] we introduced the solar foreshock coordinates developed by lGreenstadt and Bauml (|1986|) . 
and present our main results. Discussion and our conclusions are summarized in section \5\ and |6j 



2. Observations 

The observations used for the present study consist of the three components of magnetic field 
measured by Cassini MAG obtained upstream from the Kronian shock during the first fifteen 
months of Cassini's orbital data, i.e. from SOI in June 2004 through August 2005. Cassini's tra- 
jectory in Kronian Solar Orbital (KSO) coordinates for that period is shown in Figured] In the 
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KSO coordinate system centered on Saturn, the x kso axis points to the Sun, the y kso axis is anti- 
parallel to Saturn's orbital velocity and the z kso axis points towards the north pole of the Ecliptic. 
During this phase of the orbital tour, the spacecraft explored the dawn side of Saturn's magne- 
tosphere at low Kronographic latitudes. The Cassini dual magnetometer investigation consists of 
a Vector Helium Magnetometer (VHM) and a Fluxgate Magnetometer (FGM) which provide re- 
dundant, fast vector measurements of the ambient magnetic field over a wide range. The VHM 
provides accurate vector measurements with a resolution of 2 s _I over a range of +256 nT, whereas 
FGM samples the magnetic field over a larger ra nge (+65655 nT) and at higher frequency (32 s" 1 ). 
This dual technique is extensively described in iDougherty et al. For the purpose of this 

paper, we only consider the data provided by the VHM. Unfortunately, CAPS could not be used 
in this study due to the absence of periods with adequate pointing. 




Figure 1: Cassini's trajectory, since SOI in June 2004 until August 2005, projected into the x-y and the x-z plane in 
Kronocentric Solar Orbital coordinates. The average bow shock fit is the gray-solid line. 

In the absence of multi- spacecraft observations, the reconstruction of the foreshock's geometry 
requires a model of the Kronian bow shock. In this work, we used a static model bas ed on a fit of 
the crossings using a conic section. The functional form for the fit was introduced by Slavin et al. 



(|1985f) . assuming that the bow shock is axially symmetric about the solar wind flow direction. The 
general equation of a conic section is given by 

L 



1 + ecow 

where r is the distance from the planet to a point on the shock surface, 8 is the corresponding polar 
coordinate angle with respect to the symmetry axis, L is the semila tus rectum (size pa rameter) 
and e is the eccentricity. An empirical model has been presented by Ma sters et al. (2008) using 



measurements from Cassini MAG, CAPS' electron spectrometer sensor (ELS) and density mea- 
surements from the Radio Plasma Wave System (RPWS). Assuming a constant solar wind speed 
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(« 500 km/s) they found an eccentricity of e = 1 .05 + 0.09 (hyperboloid) and a size parameter 
for the average bow shock location of L = (51 + 2) R s , which implies that the average standoff 
distance is (25 ± I) Rs- Throughout their study, they did not find any dependence of the global 
configuration of the bow shock model with the IMF cone angle. 



3. Characterization of ULF waves 

During Cassini's excursions into the solar wind, we looked for intervals with low-frequency 
fluctuations in the magnetic field components. We found intervals with and without the presence 
of waves. Figure[2]and figure[3]show two examples of the types of wave events found in this paper. 

For all the wave events, we made a characterization (polarization and frequency) in the space- 
craft frame. We calculated their frequencies to (in the spacecraft frame). As a result of this study, 
we find two distinct types of oscillations with different properties, depending on whether their 
frequencies are below or above the local proton cyclotron frequency (Qh+ = eB /m H c). We 
also studied their polarization and propagat ion with respect to the ambi ent magnetic field using 



the Minimum Variance Analysis (MVA) (ISonnerup and Scheibld . [1998). The MVA consists in 



building the variance matrix in terms of the measured magnetic field components for a given 
time interval and finding the three eigenvalues and corresponding eigenvectors. The eigenvector 
corresponding to the smallest eigenvalue A3 (the maximum and intermediate eigenvalues are re- 
spectively At and A2) is used as an estimate of the direction of propagation of a plane wave. Note 
that the eigenvector set of the variance matrix provides a convenient natural coordinate system in 
which to display and analyse the data. 




Figure 2: Example of a wave train detected by the VHM on board Cassini on 17 March (day 076) 2005 between 
08:35:00 UT and 10:35:00 UT. The quasi-monochromatic packet between 08:31: 12 UT and 08:34:04 UT is analyzed 
in Figure HI 
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Figure 3: Example of a wave train detected by the VHM on board Cassini on 27 December (day 362) 2004 between 
01:00:00 UT and 03:00:00 UT. This wave event corresponds to a non linear packet. The packet between 01:30:36 UT 
and 01:33:00 UT is analyzed in Figure|5] 

Figure [2] shows an example of waves between 08:31:08 UT and 10:35:00 UT on 17 March (day 
076) 2005. We obtained that these waves are quasi-monochromatic and steepened with frequencies 
above Q H +- hi the spacecraft frame, these waves have periods of the order of 60 s. According to the 
ambient magnetic field magnitude (~0.3 nT), the period is significantly lower than the local proton 
cyclotron period T H+ ~ 300 s. Figure |4] shows the MVA results for a quasi-monochromatic wave 
packet shown in Figure [2] (interval 08:31:12-08:34:04 UT). Figured shows the components of 
the magnetic field along the maximum, intermediate and minimum variance direction. Figured 
andSt show the projection of the wave magnetic field (hodograms) on the maximum-intermediate 
and the minimum-intermediate variance planes, respectively. The circle and the asterisk indicate 
the beginning and the end of the hodogram, respectively. A large ratio between the intermediate 
and minimum eigenvalues of the variance matrix (A2M3 « 10) shows a clear minimum variance 
direction. The hodogram in |4fc shows that the polarization on the minimum variance plane is 
circular right-handed with respect to the mean magnetic field (B 3 > 0). The angle between the 
mean magnetic field and the minimum variance eigenvector is 9 kB = 20° + 1°, revealing that these 
waves propagate in a slightly oblique direction with respect to the ambient magnetic field. 

Figure |3] shows an example of a wave train seen by Cassini VHM between 01:00:00 UT and 
03:00:00 UT on 27 December (day 362) 2004. We obtained that this kind of waves are phase 
steepened with frequencies smaller than Qh+- These waves have periods of the order of 5 to 10 
min in the spacecraft frame, and they were the most frequently observed. According to the ambient 
magnetic field, 2 to 3 times the local proton cyclotron period for an IMF magnitude between 0.35 
nT and 0.5 nT. In several cases, we saw a steepening front located at the right of waves with a 
higher frequency wave packet attached to it. We note a decrease in amplitude and in the period of 
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these oscillations with increasing distance from the steepening front. Figure \5\ shows the results 
of the MVA applied on one period of these waves between 01:30:36 UT and 01:33:00 UT. From 
left to right, it can be seen how an early linear polarization is followed by a circular polarization 
toward the end of the interval. The angle between the minimum variance vector and the mean 
magnetic field suggests that the propagation of these waves is quasi-parallel (# kB = 4° + 1°) to the 
ambient magnetic field. The hodogram in Figure [5] shows that the magnetic field rotation around 
the minimum variance direction is left-handed with respect to the ambient magnetic field (B 3 <0). 
The circle and the asterisk indicate the beginning and the end of the hodogram , respe ctively. It is 
worth noticing that we identified the same two categories than iBertucci et al.l (|2007 ) had found, 
for a much bigger Cassini MAG's data set. 




Figure 4: (a) Magnetic field components along the maximum (solid line), intermediate (dashed line) and minimum 
variance direction (dot-dashed line), (b) Hodogram showing the magnetic field in the intermediate-minimum variance 
plane and (c) in the maximum-intermediate variance plane. The circle and the asterisk indicate the beginning and the 
end of the hodogram, respectively. 
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Figure 5: (a) Magnetic field components along the maximum (solid line), intermediate (dashed line) and minimum 
variance direction (dot-dashed line), (b) Hodogram showing the magnetic field in the intermediate-minimum variance 
plane and (c) in the maximum-intermediate variance plane. The circle and the asterisk indicate the beginning and the 
end of the hodogram, respectively. 



4. Determination of Saturn's ULF wave foreshock boundary: Results 

4.1. Data Selection Criterion 

To identify the ULF wave foreshock boundary using the data set described in section [21 we 
looked for beginnings (or endings) of intervals in which the magnetometer detected ULF waves. 
Due to the lack of measurements of other physical variables (solar wind density, velocity, dynamic 
pressure and composition), it is essential to establish a criterion to differentiate a ULF wave fore- 
shock boundary crossing from other possible phenomena, such as discontinuities propagating in 
the solar wind. 

For a correct determination of the ULF wave foreshock boundary, we considered those cross- 
ings in which Cassini was entering or leaving the wave region under steady IMF conditions. Using 
the whole orbital data set of Cassini's VHM, we identified a total of 59 beginnings or endings of 
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Cassini VHM - KSO Coordinates - 2005-207T03:00:00 - 2005-207T08:00:00 




Figure 6: Example of stationary crossing (entry) of the ULF wave foreshock boundary as detected by the VHM 
magnetometer on board Cassini between 03:00:00 UT and 08:00:00 UT on 26 July (day 207) 2005. Average values 
in the wave zone (Z? w ) and in the zone without waves (B nw ) are in dashed-gray line. The solid-gray lines correspond 
to (BJ>±o-p. 



intervals in which the magnetometer detected ULF waves. For this determination, we made no 
distinction regarding the kind of the field oscillations, requiring only that the transition from or 
to the wave region would be clearly apparent. We defined the following selection criterion for 
each component of the magnetic field (j = x,y,z): if the difference between the average values in 
the wave zone (BJ) and in the zone with no waves (5 nw ) is smaller than the standard deviation in 
the wave zone (crj), we consider that Cassini crossed a stationary ULF wave foreshock boundary 
(Figure©. If any of the three components did not satisfy this condition, the event discarded from 
our analysis (Figure [7]). 

The average values for each component of the magnetic field in the wave zone and in the zone 
with no waves are the mean values of the observations over two hours before (and after) each 
apparent crossing. We found that this particular time interval is sufficiently long to obtain values 
representative of the mean magnetic field, because the longest wave periods found are 10 minutes 
long. Following this selective criterion, we reduced the original data set of 59 crossings to only 
21 stationary crossings. In a second part of the study, for each component of the magnetic field 
we changed this l.aj condition, requesting the difference between BJ and 5 nw to be smaller than 
1/2. crj (more restrictive criterion) and 3/2.cr w (less restrictive criterion). 

4.2. Foreshock Coordinates 

In order to identify the ULF wave foreshock boundary independently from the changes in the 
IMF and the corresponding location of the bow shock, we analyz e our data in a particular set of 



coordinates. In the case of the Earth. iGreenstadt and Baum (1986) introduced the so-called solar 



foreshock coordinates. In this coordinate system, the x axis points toward the Sun (it is parallel to 
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Cassini VHM - KSO Coordinates - 2004-1 92T09:30:00 - 2004-1 92T1 5:30:00 




Figure 7: Example of non stationary crossing (exit) of the ULF wave foreshock boundary, as detected by the VHM 
magnetometer on board Cassini between 09:30:00 UT and 15:30:00 UT on 10 July (day 192) 2004. Average values 
in the wave zone (Z? w ) and in the zone without waves (B nw ) are in dashed-gray line. The solid-gray lines correspond 
to (BJ>±o-p. 



x kso ), and the x-y plane is the v sw -B plane which passes through the location of the spacecraft at a 
particular stationary crossing. 

For each crossing, in or out of th e ULF wave foresh ock under steady IMF conditions, we use 
the bow sho ck hyperboloidal fit f rom Sl avin et al. (1985). We used the eccentricity e = 1.05+0.09 
obtained by Masters et al. (2008) and kept it constant throughout our study. In each crossing, we 
estimate the value of L using equation (OQ) and the location of the nearest bow shock crossing. For 
this purpose, we used the list of times and locations of Cassini's bow shock crossings between 27 
June (day 027) 2004 and 12 August (day 224) 2005 published bv lMasters etaD (l2008h . The error 
in the determination of the parameter size was AL = 2 R s . 

Assuming v sw || x kso , we calculate the foreshock coordinates in the v sw -B plane or x-y plane 
(see Figure [8]) as: 

= (y - yd 

sen6> Bx 



v = 



(y a - yd 



(2) 



+ Xi 



tan# Bx 

where #b x = acos _1 (B • x kso /B) is the IMF cone angle, (jc,-, y,) is the intersection point between the 
tangent IMF line and the bow shock fit, while (x ,y ) is the observation point (Cassini's crossing 
location). It is worth mentioning that the v sw -B plane in general does not contain the x kS o axis 
(except for the particular v sw -B plane which passes through the center of the planet). In fact, 
there sh ould be a third coordinate to measure the distance of any v sw -B plane to the center of the 
planet. iGreenstadt and Baum (1986) have deliberately ignored that third coordinate, and showed 
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a remarkable correlation in the ji-v plane. As shown in Figure [H the non-orthogonal coordinate 
set (ji, v) has its origin in the tangent point (xi,yi), /u being the distance along the tangent IMF line 
and v being the horizontal distance from that particular field line to Cassini's location. 




Figure 8: Definition of the foreshock coordinates in the v sw -B plane, i.e. the x-y plane. The coordinate /i is the 
distance along the tangent magnetic field line between the tangent point and the observation point. The coordinate v 
is the distance along the Xk so direction between the tangent magnetic field line and the observation point. 

For the 21 stationary crossings considered in section I4T1 we have calculated their locations in 
terms of the coordinate set (//, v) . At Eart h, the location of the ULF wave foreshock boundary w as 



found to depend on the IMF cone angle (|Greenstadt and Bauml (I1986|) . Le and Russelll (11992ar) ). 



For this reason, we analyzed two separate sets of data, one with # Bx < 45° and the other with 
#Bx > 45°, i.e. for small and large cone angles. We performed a scatter plot of the 21 ULF wave 
foreshock boundary crossings considered. Figure[9]shows the cases with 6 Bx > 45° in black circles 
and those with 9b x < 45° in gray circles, and the straight line is our best linear fit (black line) for 
all our crossings. The small cone angle cases (gray circles) correspond to tangent lines close to 
the asymptote of the hyperbola = cos _1 (l /e) « 18°), and therefore have relatively large error 
bars. However, our best fit considering only the theta Bx > 45° cases, yields vQx) = a.fi + b 
{a = 0.47 + 0.04; b = -5.39 ± 1.10), which is indistinguishable from the result displayed in Figure 
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9. 

Saturn's ULF foreshock boundary (for all Bx cases) 




HCRs) 

Figure 9: ULF wave foreshock boundary for 0b x > 45° (dark circles) and 6b x < 45° (gray squares). If we only consider 
the 9bx > 45° cases, our best fit yields v(yu) = a./u + b (a — 0.47 + 0.04; b = -5.65 + 1.36), which is indistinguishable 
from the best fit displayed in this Figure. 

We also studied the variability of this result as our selection criterion become more stringent 
(see sub-section l4.ll) . In Figure[9]we only consider those crossings in which the three components 
of the magnetic field (J = x, y, z) satisfy that the difference between the average values in the wave 
zone and in the zone without waves is smaller than the standard deviation in the wave zone (l.crj). 
As we mentioned in sub-section EJ3 if any of the three components did not satisfy this condition, 
the crossing was discarded from our analysis. 

In a more restrictive criterion, for each component of the magnetic field, we changed the l.crj 
condition to 1/2. crj. In this way we reduced the original data set of 21 stationary crossings to only 
16 crossings. Considering only these 16 cases, our best fit yields v(jj) = a.fj. + b (a = 0.46 + 0.04; 
b = -5.26 ± 1.12), which is almost identical to the result displayed in Figure|9j We also tried with 
a less restrictive criterion, changing the original l.cr w to 3/2.cr w . In this case our original data set 
of 21 crossings, increase to 29 crossings, and our best fit yields v(p) = a.fi + b (a = 0.49 + 0.03; 
b = -5.99 + 1.01), which is again, within the errors, the same result displayed in Figure[9j 
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Saturn's ULF foreshock boundary (the three criteria) 
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Figure 10: ULF wave foreshock boundary for different criteria: the crossings in black correspond to the more re- 
strictive criterion (l/2.cr w condition); the crossings in gray are added when the l.crf criterion is considered, while the 
crossings indicated by white are added when the less restrictive criterion of 3/2. cr w is considered. 



In Figure \W\ we show the results for the three criteria. Black circles are those crossings sat- 
isfying the most stringent criterion (1/2. cr w ), gray circles correspond to crossings between 1/2. cr w 
and l.crj, while white circles are those between l.crj and 3/2. cr w . In all these cases, the uncertain- 
ties in the determination of the three magnetic field components were dominated by the statistical 
error given by the corresponding standard deviation. As mentioned, for the semilatus rectum (size 
parameter) L and the eccentricity e from the hyperboloidal model, we consider the uncertainties 
presented in section [3j 

4.3. 6 Bn = 45° Curves 

For each crossing studied, we extend the straight magnetic field line and check for connection 
to the bow shock fit. As mentioned in the p revious section, we checked for magnetic field line 



connectivity to the bow shock by using the ISlavin et al.l 41985) fit, with fixed eccentricity e = 



1.05 + 0.09. We varied L so as to make the bow shock fit to coincide with the closest bow shock 
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crossing made by Cassini. 

If we find that Cassini is connected to the shock by straight magnetic field lines, we identify 
the intersection point on the bow shock fit. Once this point has been identified, the angle between 
the magnetic field and the shock normal at the shock surface, 6 Bn , can be calculated. In all the 
wave events that we identified we found that the magnetic field line intersects the bow shock fit, 
i.e. Cassini was in the region magnetically connected to the bow shock. In particular, we find that 
the increment in wave activity is associated to foreshock regions for which 9 Bn < 45°. In fact, 12 
out of the 21 crossings have 6 Bn = 45° ± 5°, while the remaining 9 crossings have 6 Bn within the 
range between 35° and 55°. If we consider the 21 crossings, the average value is 6 Bn = 42° and the 
standard deviation is <Tg Bn = 9°. 




x/L 



Figure 11: The curves of 8 Bn = 45° (dot-dashed line) and 6 Bn = 90° (dashed line). In the v sw -B plane z = 0, the 
average bow shock fit (in units of L), the tangent field line for 6 Bx = 85° (point line), and the magnetic field line 
corresponding to 6 Bn = 45° (solid line). 

It is important to emphasize that 6 Bn = 45° constitutes a natural and conventional distinction 
between quasi-parallel and quasi perpendicular shocks. For a given magnetic field orientation and 

15 



parameter size, there are two lines on the hyperboloidal shock surface, for which 6 Bn = 45°. In 
Figure[TT] we show a a schematic view of the bow shock (in units of L) and the lines corresponding 
to Bn = 45° (dot-dashed lines) and 6 Bn = 90° (dashed lines), where the field lines are tangent to the 
hyperboloidal shock. The x-y plane in FigureQTJis parallel to the magnetic field lines and contains 
the planet in (x,y) = (0,0). Within this context, we claim that the set of magnetic field lines that 
cross the shock on the 6 Bn = 45° curves will determine the ULF wave foreshock boundary. On 
other hand, a perhaps more abstract identific ation of the ULF wave fore shock boundary is given 
by the best fit of crossings on the fi-v plane (IGreenstadt and Bauml .ll986). 




Figure 12: The (fi,v) pair for the location of points satisfying 6b„ = 45° at different (parallel) v sw -B planes for 6b x - 85° 
and parameter size L — 5 1 Rs . The distance of each plane to the planet is labeled in units of L. 

Therefore, we decided to study the consistency of our linear best fit in the fi-v plane and the 
curve corresponding to 6 Bn = 45°. Assuming a constant 6 Bx (the magnetic field lines lay in the 
x-y plane) and parameter size L = 51 R s , we computed the (ji,v) pair for the location of the point 
that satisfies 6 Bn = 45° at different (parallel) v sw -B planes (see Figure [[2]). The corresponding 
distance of each plane (in units of L) to the planet is labeled. We can clearly see in Figure [T2l an 
approximate proportionality between the /i-v coordinates at any given plane. This implies that if 
Cassini crossed right at the intersection between the bow shock and the 6 Bn = 45° field line, we 
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should expect a straight line passing through the origin in the fi-v plane (see Figure [9]). The fact 
that Cassini crosses that very field line necessarily at a distance from that intersection toward the 
upstream direction, causes a systematic increase in the value of fi, leaving v essentially unchanged. 
The distance from the bow shock measured along field lines range between 0.1 R s and 17 R s for 
the 21 crossings considered. On the one hand, this effect will reduce the quality of the linear fit, 
and on the other hand it will shift the line on the fi-v plane to the right (negative intercept), exactly 
as observed in Figure [9j 



5. Discussion 

All the wave events detected by Cassini's magnetometer show evidence of magnetic connectiv- 
ity to Saturn's bow shock, which is a clear indication that these waves are associated with Saturn's 
foreshoc k. In ab s ence o f a reliable physical model of the bow shock, we used the empirical model 



of Masters et al.1 (|2008|) . which is based on extensive observations from Cassini made during the 
same time span of our work (15 months from June 2004 through August 2005). Assuming that 
the bow shock is axially symmetric about the solar wind flow direction, the average bow shock 
location is obtained by fitting a conic section to the bow shock crossings using a nonlinear least 
squares technique. This fit provides an adequate description of the average shape, given by an ec- 
centricity of e = 1.05+0.09, and also that of its average location, corresponding to L = (5l±2)R s . 
In our study, we assume that the shape of Saturn's bow shock remains constant throughout our ob- 
servational time span, while its location varies in response to time variations of the solar wind's 
ram pressure. In other words, for each wave event or ULF wave foreshock boundary crossing, we 
assumed e to remain constant. Meanwhile, the semilatus rectum L, which determines the distance 
of the bow shock to the planet, was allowed to change to account for this pressure effect. We 
found values of L between 53 + 2 R s and 78 + 2 R s . We also calculated the parameter size using 
an average between the two closest bow shock crossing locations for a given event (wave or cross- 
ing). For the sake of comparison, we also made our calculation assuming that the parameter size 
remains constant (L = 62 R s , which is the average value) for all the wave events. We found, that 
the bow shock fit was not representative of many of the upstream events studied, since they fell in 
the downstream region of this static fit. Therefore, we concluded that the bow shock fit calculated 
using equation ([T]) with a varying L parameter gives a more accurate empirical representation of 
the real Kronian bow shock. 

For all the wave events observed, we performed a preliminary characterization based on their 
frequency and polarization as measured in the spacecraft reference frame. According to their 
frequencies, we observe two different types of oscillations. Waves with frequencies below Q H+ 
are the most frequently observed, they are phase steepened with periods of the order of 5 to 10 
minutes in the spacecraft frame. As we mentioned in section|3j we saw steepening fronts located at 
the right of waves with a higher frequency wave packets attached to them. We also note a decrease 
in amplitude and in the period of the oscillations within each wave packet with increasing distance 
from the steepening front. The angle between the minimum variance vector and the mean magnetic 
field suggests that the propagation of these waves is quasi-parallel. The magnetic field rotation 
around the minimum variance di rection is left-handed with respect to the IMF direction. As was 



discussed in Be rtucci et al.l (|2007|) . these signatures suggest that these waves are ion/ion resonant 
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right-hand (fast magnetosonic) mode waves which steepen during the nonlinear regime and emit 
a dispersive whistler to stop the steepening. We also observe waves with frequencies above Q.u+, 
which are either quasi-monochromatic or steepened, with periods of the order of 60 s (significantly 
lower than the local proton cyclotron period T H+ ~ 300 s). For the quasi-monochromatic events, 
their corresponding hodograms show that the polarization is circular right-handed with respect 
to the mean magnetic field. The angle between the ambient magnetic field and the minimum 
variance direction reveals that these waves propagate in a slightly oblique direction with respect 
to the IMF orientation. This kind of wave packets is an exception rather than the rule (we found 
ten events only). At the Earth's foreshock, there is a partic ular type of ULF wave s, the so-called 
30 s (period) waves which are quasi-monochromatic (e.g. iLe and RussellL Il992bh . This kind of 
waves are found always near to the ULF wave foreshock boundary. In our study we did not find 
quasi-monochromatic waves next to the Saturnian ULF wave foreshock boundary. A more detailed 
statistical study would shed some light about the nature of this second group of waves and whether 
they are the Saturnian equivalent of the 30 s modes found in the Earth's foreshock. 

When Cassini is magnetically connected to the bow shock fit by straight magnetic field lines 
(since we are assuming a uniform external magnetic field), we identify the intersection point be- 
tween the field line and the bow shock fit as well as the angle 6 Bn . In all the crossings that we 
observed, we find that the presence of magnetic fluctuations is associated with values of 9 Bn cor- 
responding t o quasi-parallel geome tries (0g n < 45°). This result is fully consistent with the one 
obtained by Le and Russelll (|1992a|) at Earth. Using ISEE data, they examined the presence of 
ULF waves in the region immediately upstream from the Earth's bow shock. Their statistical study 
shows a well-defined boundary on the v sw -B plane, separating the disturbed from the undisturbed 
magnetic field at 6 Bn ~ 50°. Our result for Saturn's foreshock is consistent with this scenario, 
regardless of whe ther these waves are immediately adjacent to the bow shock or deep into the 
foreshock region. ISchwartz et al. ( 1983|) developed a theoretical framework for studying the gen- 
eral trajectories of ions reflected or leaked upstream from the Earth's bow shock and subjected 
only to the Lorenz force in a steady IMF. In the case known as specular reflection, the ions bounce 
off the shock potential reversing their velocity component along th e shock normal h, while main - 
taining their parallel component. In the de Hoffman-Teller frame (Ide Hoffman and Te ller. 1950), 
assuming that the bow shock is infinitesimally thin and has no overshoot. ISchwartz et al] (|1983|) 
found that after interacting with the shock, the guiding center motion of incident solar wind ions 
is reflected to the upstream region with 6 Bn < 45°. For 6 Bn > 45°, specularly reflected ions always 
re-encounter the bow shock with sufficient energy to penetrate into the downstream region. It is 
worth mentioning that these specularly reflected ions oriented to the upstream region for quasi- 
parallel geometries would have no conserving magnetic moment, in c ontrast with the field-alig ned 
reflected beams, for which there is a magnetic moment conservation (IPaschmann et all Il980l) . In 
the absence of Cassini plasma observations that could confirm the presence of backstreaming ions 
giving rise to an ion foreshock region, we cannot observationally associate the presence of waves 
with specific ion distributions. However, if the mechanism of generation of these waves is local, 
i.e. the wave growth rate y is sufficiently large, we could expect specularly reflected ions to be 
responsible for these waves since they are the only population whose spatial distribution depends 
on 9 Bn . In this case, for quasi-parallel geometries, the reflected ions can go upstream from the 
bow shock and be detectes through ULF wave activity. Nevertheless, a full discussion on the wave 

18 



generation caused by particle distributions as a function of 6 Bn , is beyond the purpose of this study. 

For a correct determination of Saturn's ULF wave foreshock boundary, we only considered 
crossings in which Cassini was entering or leaving the wave zone under steady IMF conditions, 
i.e. stationary crossings. We defined a criterion based on the variability of the magnetic field in 
order to make a proper identification of these stationary crossings. As a result, we reduced our 
original list of crossings to o nly 21 stationary ones. At Earth, using single spacecraft techniques, 



Greenstadt and Baum ( 1986) obtained different orientations of ULF wave foreshock boundary for 



different cone angles (more specifically, they considered two ran ges: 6r x = 20° - 30° an d 6 Bx = 40° 
- 50°). In agreement with Parker's angle prediction for Saturn ( Jackman et all |2008|) . our study 
shows that 83% of the crossings correspond to cone angles larger than 55°. Therefore, the range 
of 9 Bx considered is different to that at Earth. Nevertheless, we analyzed the correlation between 
ii-v coordinates for two ranges of 6 Bx : one set of data with 6 Bx < 45° and the other with 9 Bx > 45°. 
A comparison of the best linear fit for the two data set indicates no significant difference between 
them (see Figure [9]). 

This result, however, does not rule out a dependence on 6 Bx since the errors in the determination 
of fi-v become larger with decreasing Bx . More specifically, the main source of error propagation 
in the determination of each {p., v) pair is the intersection point between the tangent IMF line and 
the bow shock fit. This intersection point depends critically on the IMF cone angle, which is 
therefore the main source of error. In particular, small cone angles correspond to almost tangent 
magnetic field lines close to the asymptotes of the hyperbola. On the other hand, the appearance 
of waves is observed to occur at 6 Bn values of ~ 45°. The surface determined by the condition 
6 Bn = 45° is mapped on the ii-v plane as straight line, as illustrated by Figure [121 As a result, 
we find that there is an equivalence between the ULF wave foreshock boundary determined by 
6 Bn ~ 45° and the best fit in the fx-v plane citepGB1986. They are in fact two alternative ways to 
determine the same spatial region within the foreshock. 

Finally, we studied the robustness of the ii-v correlation as our selection criterion becomes 
more stringent. For instance, we change the I. erf condition to a more stringent value of l/2.crf 
for the three components of magnetic field, or change it into the less restrictive value of 3/2. cr w for 
each component. We find that regardless of whether we choose a more or less restrictive criterion, 
our best fit remains the same, within the fit errors. We also find that for these three different 
criteria, the boundary does not seem to depend on whether the cone angle is large or small. 



6. Conclusions 

Using Cassini's data from SOI in June 2004 through August 2005, we conducted a detailed 
survey and analysis of ULF waves upstream from Saturn's bow shock. All the wave events show 
evidence of magnetic connectivity to Saturn's bow shock, therefore we conclude that these wave 
events are undoubtedly associated with Saturn's foreshock. 

As for their frequencies, we identify two distinct types of wave populations. The most fre- 
quently observed have frequencies below Q H + and are phase steepened, with periods of the order of 
5 to 10 minutes in the spacecraft frame. These waves are often accompanied by precursor whistler 
wave trains. In agreement with previous works, we suggest that these waves are ion/ion reso- 
nant right-hand (fast magnetosonic) mode waves which steepen during the nonlinear regime and 
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emit a dispersive whistler to stop the steepening. We also observe waves with frequencies above 
Q.u+, which appear as either quasi-monochromatic or steepened, with periods of ~ 1 minute. The 
quasi-monochromatic events have a circular right-handed polarization with respect to the mean 
magnetic field and a propagation slightly oblique with respect to the IMF. A more detailed study 
will be needed in order to conclude whether the latter are the Saturnian equivalent of the 30 s 
modes found at Eart h. It is worth noticing that we have identified the same two categories than 
Bertucci et al.1 (12007|) have found, although for a much bigger Cassini MAG's data set. 

We identified 21 stationary crossings inbound and outbound from the ULF wave foreshock 
region. We calculated their solar foreshock coordinates in the v sw -B plane and we have identified 
for the first time Saturn's ULF wave foreshock boundary. In the fi-v plane we do not find a clear 
dependence between the foreshock boundary and the IMF cone angle. We also found that the 
presence of waves is associated with the change in 6 Bn to quasi-parallel geometries. Moreover, we 
find that our determination of the ULF wave foreshock boundary as the surfac e given by 6r„ = 45°, 
is ind eed consistent with the linear fit in the fi-v plane, first proposed by iGreenstadt and Baum 
dl986h for the Earth's ULF wave foreshock. In this regard, we speculate that the specular reflection 
is the candidate process for the reflected ions (gyrophase-bunched and diffuse distributions), since 
it is the only process to produce a signature in Finally, we studied the robustness of the fi-v 
correlation for a more restrictive criterion. We found that regardless of the criterion that we use, 
we obtain the same ULF wave foreshock boundary in the fi-v plane. 



References 

Achilleos, N., Bertucci, C, Russell, C. T., Hospodarsky, G. B., Rymer, A. M., Arridge, C. S., Burton, M. E., 
Dougherty, M. K., Hendricks, S., Smith, E. J., Tsurutani, B. T., 2006. Orientation, location, and velocity of Saturn's 
bow shock: Initial results from the Cassini spacecraft. J. Geophys. Res. Ill, A03201. 

Acuna, M. H., Ness, N. E, 1980. The magnetic field of Saturn: Pioneer 11 observations. Science 207, 444-446. 

Balogh, A., Dunlop, M. W., Cowley, S. W. H., Southwood, D. J., Thomlinson, J. G, Glassmeier, K. H., Musmann, 
G, Liihr, H., Buchert, S., Acuna, M. H., Fairfield, D. H., Slavin, J. A., Riedler, W., Schwingenschuh, K., Kivelson, 
M. G, 1997. The Cluster magnetic field investigation. Space Sci. Rev. 79, 65-91. 

Bavassano-Cattaneo, M. B., Moreno G, Scotto, M. T, Acuna, M., 1987. Observations of Large-Amplitude MHD 
Waves in Jupiter's Foreshock in Connection With a Quasi-Perpendicular Shock Structure. J. Geophys. Res. 
92(A11), 12413-12418. 

Bavassano-Cattaneo, M. B., Cattaneo, P., Moreno, G, Lepping, R. P., 1991. Upstream waves in Saturn's foreshock. 

Geophys. Res. Lett. 18, 797-800. 
Belcher, J. W., Richardson, J. D., Vasyliunas, V. M., Lepping, R. P., Ness, N. E, Smith, C. W., 1991. Low-frequency 

waves in the solar wind near Neptune. Geophys. Res. Lett. 18, 1071-1074 
Bertucci, C, Achilleos, N., Mazelle, C, Hospodarsky, G. B., Thomsen, M., Dougherty, M. K., Kurth, W., 2007. 

Low-frequency waves in the foreshock of Saturn: First results from Cassini. J. Geophys. Res. 1 12, A09219. 
Blomberg, L., Cumnock, J., Glassmeier, K.-H, Treumann, R., 2007. Plasma Waves in the Hermean Magnetosphere. 

Space Sci. Rev. 132, 575-591. 
Convery, P. D., Gary, S. P., 1997. Electromagnetic proton cyclotron ring instability: Threshold and saturation. J. 

Geophys. Res. 102, 2351-2358. 
de Hoffman, E, Teller, E., 1950. Magneto-hydrodynamic shocks. Phys. Rev. 80, 692-703. 

Diodato, L., Greenstadt, E. W., Moreno, G, Formisano, V., 1976. A Statistical Study of the Upstream Wave Boundary 

Outside the Earth's Bow Shock. J. Geophys. Res. 81(1), 199-204. 
Dougherty, M. K., Kellock, S., Southwood, D. J., Balogh, A., Smith, E. J., Tsurutani, B. T, Gerlach, B., Glassmeier, 

K. H, Gleim, E, Russell, C. T, Erdos, G, Neubauer, F. M., Cowley, S. W. H, 2004. The Cassini Magnetic Field 

Investigation. Space Sci. Rev., 114, 331-383. 

20 



Gary, S. P., 1993. Theory of Plasma Microestabilities, Cambridge Univ. Press, New York. 

Eastwood, J. P., Balogh, A., Lucek, E. A., Mazelle, C, Dandouras, I., 2003. On the existence of Alfven waves in the 

terrestrial foreshock. Ann. Geo. 21, 1457-1465. 
Eastwood, J. P., Lucek, E. A., Mazelle, C, Meziane, K., Narita, Y., Pickett, J. and Treumann, R. A., 2005. The 

Foreshock. Space Sci. Rev. 118, 41-94. 
Fairfield, D. H., K. W. Behannon, 1976. Bow Shock and Magnetosheath Waves at Mercury. J. Geophys. Res. 81(22), 

3897-3906. 

Fuselier, S. A., Gurnett, D. A., Fitzenreiter, R. J., 1985. The downshift of electron plasma oscillations in the electron 
foreshock region. J. Geophys. Res. 90, 3935. 

Greenstadt, E., Baum, L., 1986. Earth's Compressional Foreshock Boundary Revisited: Observations by the ISEE 1 
Magnetometer. J. Geophys. Res. 91(A8), 9001-9006. 

Greenstadt, E. W., Baum, L. W., Jordan, K. F, Russell, C. T, 1987. The Compressional ULF Foreshock Boundary of 
Venus: Observations by the PVO Magnetometer. J. Geophys. Res. 92(A4), 3380-3384. 

Gosling, J. T, Asbridge, J. R., Bame, S. J., Paschmann, G., Sckopke, N., 1978. Observations of two distinct popula- 
tions of bow shock ions in the upstream solar wind. Geophys. Res. Lett. 5, 957-960. 

Gurgiolo, C, Parks, G. K., Mauk, B. H., 1983. Upstream gyrophase bunched ions: a mechanism for creation at 
the bow shock and the growth of velocity space structure through gyrophase mixing. J. Geophys. Res. 88(11), 
9093-9100. 

Gosling, J. T, Thomsen, M. F, Bame, S. J., Feldman, W. C, Paschmann, G., Sckopke, N., 1982. Evidence for 

specularly reflected ions upstream from the quasi-parallel bow shock. Geophys. Res. Lett. 9(12), 1333-1336. 
Hoppe, M. M., Russell, C. T, 1983. Plasma rest frame frequencies and polarizations of the low-frequency upstream 

waves: ISEE 1 and 2 observations. J. Geophys. Res. 88, 2021-2028. 
Hoppe, M. M., Russell, C. T, Frank, L. A., Eastman, T. E., Greenstadt, E. W., 1981. Upstream hydromagnetic waves 

and their association with backstreaming ion populations: ISEE 1 and 2 observations. J. Geophys. Res. 86, 4471. 
Hoshino, M., Terasawa, T, 1985. Numerical study of the upstream wave excitation mechanism. Nonlinear phase 

bunching of beam ions, J. Geophys. Res., 90(1), 57-64. 
Jackman, C. M., Forsyth, R. J., Dougherty, M. K., 2008. The overall configuration of the interplanetary magnetic field 

upstream of Saturn as revealed by Cassini observations. J. Geophys. Res. 113, A081 14. 
Krimigis, S. M., Sards, E. T, Mitchell, D. G., Hamilton, D. C, Bertucci, C, Dougherty, M., 2005. Spatial distribution, 

composition and charge state of energetic ions upstream from the Kronian magnetosphere. Eos Trans. AGU 86(52), 

Fall Meet. Suppl. Abstract SH53B-08. 
Le, G., Russell, C. T, 1992a. A study of ULF wave foreshock morphology - 1: ULF foreshock boundary. Planet. Space 

Sci.. 40(9), 1203-1213. 

Le, G., Russell, C. T, 1992b. A study of ULF wave foreshock morphology - II: spatial variation of ULF waves. Planet. 

Space Sci. 40(9), 1215-1225. 
Masters, A., Achilleos, N., Dougherty, M. K., Slavin, J. A., Hospodarsky, G. B., Arridge, C. S., Coates, A. J., 2008. 

An empirical model of Saturn's bow shock: Cassini observations of shock location and shape. J. of Geophys. Res. 

113,A10210. 

Mazelle, C, Le Queau, D., Meziane K., 2000. Nonlinear wave-particle interaction upstream from the Earth's bow 
shock. Nonlinear Proc. Geophys. 7, 185-190. 

Mazelle, C, Meziane, K., Le Queau, D., Wilber, M., Eastwood, J. P., Reme, H., Sauvaud, J. A., Bosqued, J. M., 
Dandouras, I., McCarthy, M., Kistler, L. M., Klecker, B., Korth, A., Bavassano-Cattaneo, M. B., Pallocchia, G., 
Lundin, R. and Balogh, A. (2003), Production of gyrating ions from nonlinear wave-particle interaction upstream 
from the Earth's bow shock: A case study from Cluster-CIS, Planet. Space Sci., 51, 785-795 

Mazelle, C, Winterhalter, D., Sauer, K., Trotignon, J. G, Acuna, M. H, Baumgartel, K., Bertucci, C, Brain, D. A., 
Brecht, S. H., Delva, M., Dubinin, E., Oieroset, M., Slavin, J., 2004. Bow shock and upstream phenomena at Mars. 
Space Sci. Rev. Ill, 115-181. 

Meziane, K., C. Mazelle, Lin, R. P., Le Queau, Larson, D. E., Parks, G. K., Lepping, R. P., 2001. Three-dimensional 
observations of gyrating ion distributions far upstream from the Earth's bow shock and their association with low- 
frequency waves. J. Geophys. Res. 106(4), 5731-5742. 

Ogilvie, K. W., von Rosenvinge, T, Durney, A. C, 1977. International Sun-Earth Explorer: A Three-Spacecraft 



21 



Program. Science 198(4313), 131-138. 
Paschmann, G., Sckopke, N., Bame, S. J., Asbridge, J. R., Russell, C. T., Greenstadt, E. W., 1979. Association of 

low-frequency waves with suprathermal ions in the upstream solar wind. Geophys. Res. Lett. 6, 209-212. 
Paschmann, G, Sckopke, N., Asbridge, J. R., Bame, S. J., Gosling, J. T., 1980. Energization of Solar Wind Ions by 

Reflection From the Earth's Bow Shock. J. Geophys. Res. 85(A9), 4689-4693. 
Paschmann, G, Sckopke, N., Papamastorakis, I., Asbridge, J. R., Bame, S. J., Gosling, J. T., 1981. Characteristics of 

reflected and diffuse ions upstream from the Earth's bow shock. J. Geophys. Res., 86, 4355-4364. 
Russell, C. T., Lepping, R. P., Smith, C. W., 1990. Upstream Waves at Uranus. J. Geophys. Res. 95(A3), 2273-2279. 
Schwartz, S. J., Thomsen, M. E, Gosling, J. T., 1983. Ions Upstream of the Earth's Bow Shock: A Theoretical 

Comparison of Alternative Source Populations. J. Geophys. Res 88(A3), 2039-2047. 
Sigsbee, K., Kletzing, C. A., Gurnett, D. A., Pickett, J. S., Balogh, A., Lucek, E., 2004. Statistical behavior of 

foreshock Langmuir waves observed by the Cluster Wideband Data Plasma Wave Receiver. Ann. Geophys. 22, 

2337-2344. 

Sigsbee, K., Kletzing, C. A., Gurnett, D. A., Pickett, J. S., Balogh, A., Lucek, E., 2004. The dependence of Langmuir 

wave amplitudes on position in Earth's foreshock. Geophys. Res. Lett. 31, L07805. 
Slavin, J. A., Smith, E. J., Spreiter, J. R., Stahara, S. S., 1985. Solar wind flow about the outer planets: Gas dynamic 

modeling of Jupiter and Saturn bow shock. J. Geophys. Res. 90, 6275. 
Sonnerup, B. U. O., Scheible, M., 1998. Minimum and Maximum Variance Analysis, in: Paschmann, G., Daly, P. W. 

(Eds), Analysis Methods for Multi-Spacecraft Data. Int. Space Sci. Inst., Bern, Switzerland, pp. 185-220. 
Thomas, B. T, Smith, E. J., 1980. The Parker spiral configuration of the interplanetary magnetic field between 1 and 

8.5 AU. J. Geophys. Res. 85, 6861-6867. 
Thomsen, M. E, DiLorenzo, J. P., McComas, D. J., Young, D. T, Crary, E J., Delapp, D., Reisenfeld, D.B., Andre, 

N., 2007. Assessment of the magnetospheric contribution to the suprathermal ions in Saturn's foreshock region. J. 

Geophys. Res. 112, A05220. 

Trotignon, J. G., Skalsky, A. ,Grard, R., Nairn, C, Klimov, S., 1992. Electron Density in the Martian Foreshock as a 
By-Product of the Electron Plasma Oscillation Observations. J. Geophys. Res. 97(A7), 10831-10840. 

Tsurutani, B. T, Rodriguez, P., 1981. Upstream waves and particles - An overview of ISEE results. J. Geophys. Res. 
86,4317. 

Went, D. R., Hospodarsky, G. B., Masters, A., Hansen, K. C, Dougherty, M. K., 2011. A new semiempirical model 
of Saturn's bow shock based on propagated solar wind parameters. J. Geophys. Res. 116, A07202. 



22 



